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Abstract 

A new dynamical model for the propagation of strangelets through the ter- 
restrial atmosphere is proposed. 

Strangelets are small lumps of strange quark matter (SQM) which consist of 
roughly equal numbers of up, down and strange quarks. Since the suggestion of 
WittenJTJjin 1984 that these strangelets, and not ordinary nuclear matter, represent 
the true ground state of QCD, this area of research has expanded considerably. The 
existence of stable or metastable lumps of SQM would have numerous implications 
for physics and astrophysics. Most importantly, they can account for the cosmolog- 
ical dark matter problem to a large extent, if not entirely 0. Moreover, the mere 
existence of strangelets may be one of the most perfect evidences for the QGP phase 
transition that is believed to have occurred in the early universe or inside neutron 
stars. 
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To establish the existence of strangelets, it is, of course, necessary to detect them. 
The obvious place to search for them would be in the cosmic rays. To this end, it 
is necessary to understand the propagation of strangelets through the earth's atmo- 
sphere, as these would have to traverse the whole atmospheric depth before arriving 
at the detectors. In this letter, we attempt to provide a realistic dynamical model 
which describes the propagation of strangelets. 

The fact that SQM is absolutely stable does not contradict the ordinary experience 
that matter is, for the most part, made of ordinary nuclear matter. This is because, 
in order for ordinary matter to change to SQM, a large number of u and d quarks 
need to change to s quarks, which would require a very high order weak interaction 
and therefore can be considered to be highly improbable [|l|. On the other hand, the 
existence of ordinary nuclear matter shows that quark matter consisting of only u and 
d quarks is unstable Q. Introduction of a third flavour reduces the energy relative 
to a two-flavour system since another Fermi well is now present and this makes the 
system stable. In the earlier works pi , calculations along the lines of Fermi gas model 
seemed to suggest that although strangelets with high atomic numbers would be 
absolutely stable, the same could not be said for strangelets with low A; the surface 
effect which must be taken into account for small strangelets would make the energy 
per baryon increase with decreasing baryon number, turning the SQM unstable. In 
later calculations [|], |J, the Fermi gas model was replaced by the spherical MIT 
hadronic bag models and the stability of strangelets with low A was examined in 
the light of this model. The works done by various authors || [7j tell us that for 
low s-quark masses, shell-like structures occur for A = 6, 18, 24, 42, 54, 60, 84, 102 
etc. These values may change somewhat with the change in the strangeness fraction 

Strangelets are expected to possess a small positive electric charge. They would 
have been neutral, if the ground state composition consisted of equal numbers of 
quarks of the three types of quark flavours, which is the most favourable state. Actu- 
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ally, however, the s quark is heavier than the other two flavours and so their number 
is slightly less than that of the u or d quarks. Hence the fortuitous cancellation 
\ n u — \ n d — \n s does not occur exactly, as a result of which we are left with a small 
residual positive charge. It can also be argued from mere experience that they cannot 
be negatively charged because although overall charge neutrality can be maintained 
by covering them up by a cloud of positrons, they will eat up every piece of ordinary 
matter they find in their path since it is energetically favorable for nuclear matter to 
convert to SQM. (For a recent survey of the physics of SQM, see f|.) On the other 
hand, a small positive charge would be helpful, as in that case, ordinary nuclei will 
be electrostatically repelled when the strangelets are moving slowly; when they move 
fast enough (highly energetic), nothing can prevent them from absorbing neutrons 
and becoming more and more tightly bound. 

It is generally assumed that the strangelets which come to the upper layer of the 
atmosphere have baryon number A ~ 1000 or more. This assumption is usually made 
on the basis of the following experimental and theoretical considerations. Some of 
the reports which suggest direct candidates for SQM are tabulated below. 



Event 


Mass 


Charge 


Counter experiment [[J 


A ~ 350 - 450 


14 


Exotic Track JK| 


A ~ 460 


20 


Price's Event 11 


A > 1000 


46 


Balloon Experiments |12, |i~3|j 


A ~ 370 


14 



The theoretical tip is usually taken from the generally accepted mass number for 
strangelets which is given in the above paragraph. It is often remarked that that 
the observation of these candidates at such large atmospheric depths (~ 500 g/cm 2 ) 
requires unusual penetrability of these lumps which means that their cross sections 
should be very small and hence the geometric size much smaller than typical nuclear 
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size. This, however, is not borne out by the models of SQM where the density of 
SQM is seen to be not much larger than ordinary nuclear matter M. 



Wilk et al. |14| , [L5| [L6|, [L7| proposed a mechanism by which the strangelets are 
able to cover great atmospheric depths. They assumed that the mass and hence 
the cross sections of the strangelets incident on the upper layers of the atmosphere 
with initial masses of the order of 10 3 a.m.u. decrease rapidly due to their collisions 
with air molecules in which a mass equal to that of the nucleus of an atmospheric 
atom is ripped off from the strangelet in every encounter. They also proposed that 
there should exist a critical mass m cr i t such that when the mass of the strangelet 
evolving out of an initially large strangelet drops below that, it simply evaporates 
into neutrons and that is what would fix the lower limit of the altitude upto which a 
strangelet would be able to penetrate. 

Let us recapitulate the basic conclusions that can be derived from the earlier works. 
Firstly, strangelets observed at the mountain altitudes typically have mass around 
300 to 450 and charge between 10 to 20. But the experimental results obtained till 
date are inconclusive and hence they do not impose a strict bound on the mass and 
charge of strangelets that can be observed in future experiments. Secondly, although 
the correlation between penetrability and geometric cross sections is usually valid for 
ordinary nuclei, the same cannot be easily extrapolated to the case of strangelets since 
these are very tightly bound massive objects and are not expected to break up as a 
result of nuclear collisions. Indeed, in a typical interaction between a strangelet and 
the nucleus of an atmospheric atom, it is more probable for the strangelet to absorb 
neutrons so that the colliding nucleus, and not the strangelet, is likely to break up 



most of the time. Hence the scheme proposed in | I4]| , namely that the mass of a 



strangelet decreases in every encounter, seems to be unrealistic. We, in this letter, 
propose an alternative scheme based on the following premises : 

1. The collision of a lump of SQM with ordinary matter results in the absorption 
of the neutrons from the colliding nucleus, as a result of which the mass of the 
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strangelet increases in every collision and it becomes more tightly bound. 

2. The initial masses of the strangelets are assumed to be small in order to obtain 
final baryon numbers which are nearly equal to the observed ones at mountain 
altitudes. The discussion above indicates that it is quite possible to have stable 
lumps of SQM with low mass numbers. This would also facilitate a somewhat 
larger flux in the cosmic rays. 

3. The speed, and hence the kinetic energy of these particles, must be such that 
they would arrive at a distance of 25 km above the sea level, surmounting 
the geomagnetic effects. We start with such an altitude since the atmospheric 
density above 25 km is low enough to be neglected. The charge of the strangelet 
is also fixed by this assumption, corresponding to a certain strangeness fraction. 

The simple assumptions proposed above give a picture more or less in accord with the 
observation of the propagation of the strangelets in the atmosphere, which can give 
useful indications of the type of things to be expected in an actual experiment. The 
description of the model is given next. We consider a situation in which a strangelet 
with a low baryon number enters the upper layers (~ 25Km from the sea level) of 
the atmosphere. To arrive at this point, a charged particle must possess a speed 
determined by the formula (see, e.g, []18| ). 

pc M cos A $ , . 

Ze r o (yi + cosh + i) 

where Mis the magnetic dipole moment of the Earth, r Q the radius of the Earth 
and is the (geomagnetic) latitude of the point of observation (~ 30°, which might 
represent a location in north eastern India), p and Ze represent the momentum and 
charge, respectively, of the particle. The magnetic field of the Earth is taken to be 
equivalent to that due to a magnetic dipole of moment M — 8.1 x 10 22 J/T, located 
near the centre of the earth, the dipole axis pointing south. We have fixed the mass, 
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initial speed and charge to be 64 amu 6.6 x 10 7 m per sec and 2 (electron charge), 
respectively, at the initial altitude of 25 km, for the purpose of the present work. 

In the course of its journey, the strangelet comes in contact with air molecules, 
mainly N 2 . During such collisions, the strangelet absorbs neutrons from some of these 
molecules, as a result of which it becomes more massive. The effect of such encounters 
is summarized in the formula 

dm s _ f x m N 

dh ~ A 1 ] 

where mj is the mass of the strangelet, the total mass of the neutrons in the 
atmospheric atom, A the mean free path of the strangelet in the atmosphere and h 
the path length traversed. (It should be emphasized here that the strangelets would 
preferentially absorb neutrons, as protons would be coulomb repelled. Nonetheless, 
some protons could still be absorbed initially when the relative velocity between the 
strangelet and the air molecule is large, leading to increase in both the mass and the 
charge of the strangelet. We do not address this issue in the present work for the 
sake of simplicity, although it can be readily seen that the rate of increase in mass 
would obviously be faster than that in charge.) In the above equation, A depends 
both on h, which determines the density of air molecules and the instantaneous mass 
of the strangelet, which relates to the interaction cross section. The mean free path 
decreases as lower altitudes are reached since the atmosphere becomes more dense 
and the collision frequency increases. Finally the factor / determines the fraction 
of neutrons that are actually absorbed out of the total number of neutrons in the 
colliding nucleus. The expression for this factor has been determined by geometric 



considerations |19j and is given by 

V 



/= | (1 -,)Wlzii) , _>[3(l-,)V»_ 1 ]fiz£V (3) 



In eqn(3), /i = Ri h +R2 and v = R -^ R - where b is the impact parameter. R 1 and R 2 are 
the radii of the strangelet and the nucleus of the atmospheric atom, respectively. / is 
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initially small but grows larger and reaches the limiting value 1 when the strangelet 
grows more massive. 

The above considerations lead us to a set of differential equations of the form 

— = -g + (v x B) -f- 4 

at m-s m s d* 

In eqn(4), — g represents the acceleration due to gravity, B is the terrestrial magnetic 
field, q = Ze and v represents the velocity of the strangelet. These equations were 
solved by the A th order Runge Kutta Method for the set of initial conditions described 
above. 

The results are shown in figs. 1-5. Figure 1 shows the variation of Altitude 
with time, the zero of time being at 25 km. The time required to reach a place which 
is about 3.6 km above the sea level (height of a typical north east Indian peak like 
'Sandakphu', where an experiment to detect strangelets in cosmic rays using a large 
detector array is being set up pOp is indicated in the figure. The next figure (Fig. 2) 



shows the variation of the mass of the strangelet with time and figure 3 shows the 
variation of the strangelet mass with altitude. It can be seen from the figures that 
the expected mass at the aforementioned altitude comes out to be about 340 amu or 
so. Figure 4 shows the variation of the mean free time with altitude. Finally, figure 5 
shows the variation of (3 = v/c with time, showing the expected decrease of the speed 
with time, which also justifies the adequacy of the nonrelativistic treatment that we 
have applied. 

In this letter, we have proposed a dynamical model of the propagation of strangelets 
through the atmosphere of the Earth. Although the model is based on simple assump- 
tions, it is realistic enough to include the difference in the interaction process which 
is expected when SQM, and not ordinary nuclei, collide with the atmospheric nu- 
clei. The effects of Earth's gravitational and magnetic fields are also included in the 
equations of motion so that meaningful information can be extracted directly from 
the resulting trajectory. The main conclusion of the model is that the exotic cosmic 
ray events with very small Z/A ratios at mountain altitudes could result from SQM 
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Figure 1: Variation of altitude with time. The arrow corresponds to an altitude of 
3.6 km from the sea level. 
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Figure 2: Variation of mass with time. The time indicated by the arrow is the time 
taken to reach an altitude 3.6 km from the sea level, starting from 25 km. 
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Figure 3: Variation of mass with altitude. The arrow corresponds to an altitude of 
3.6 km from the sea level. 
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Figure 4: Change of mean free time with altitude. 
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Figure 5: Change of (3 — v/c with time. 
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droplets (strangelets) which need not be too large initially. Thus the flux of strangelets 
in the cosmic rays may indeed be appreciable enough to make their detection by a 
large area detector at mountain altitudes a real possibility. 

The works of SB and SKG were supported in part by the Council of Scientific & 
Industrial Research, Govt, of India, New Delhi. 
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